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RBC and RBC/UKQCD Collaborations 

We report RBC and RBC/UKQCD lattice QCD numerical calculations of nucleon electroweak 
matrix elements with dynamical domain-wall fermions (DWF) quarks. The first, RBC, set of 
dynamical DWF ensembles employs two degenerate flavors of DWF quarks and the DB W2 gauge 
action. Three sea quark mass values of 0.04, 0.03 and 0.02 in lattice units are used with about 200 
gauge configurations each. The lattice cutoff is a -1 ~ 1.7GeV and the spatial volume is about 
(1.9fm) 3 . Despite the small volume, the ratio of the isovector vector and axial charges g A /g v 
and that of structure function moments (x}„_^/ (x)A«-Ad are m agreement with experiment, and 
show only very mild quark mass dependence. The second, RBC/UK, set of ensembles employs 
one strange and two degenerate (up and down) dynamical DWF quarks and Iwasaki gauge action. 
The strange quark mass is set at 0.04, and three up/down mass values of 0.03, 0.02 and 0.01 
in lattice units are used. The lattice cutoff is a ~ 1.6GeV and the spatial volume is about 
(3.0fm) 3 . Even with preliminary statistics of 25-30 gauge configurations, the ratios g A /g v and 
{ x )u-d I ( x )au-Ai1 are consistent with experiment and show only very mild quark mass dependence. 
Another structure function moment, d\, though yet to be renormalized, appears small in both sets. 
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1. Nucleon structure on the lattice 

Nucleon isovector vector and axial charges, g v and g A , are defined in neutron j3 decay form 
factors: g v = lim q2 ^ g v (q 2 ) with 

(n\V-(x)\p) = iu n [Yiigv(<l 2 )+<lxGXLi8M(<j Z )] u P e ~ m > 
and g A = lim CJ 2^ g A (q 2 ) with 

(n\A-(x)\p) = iunjs^g^+q^g^q^upe-^. 

Their ratio, g A /g v , is very accurately measured as 1.2695(29) [JTJ| - 

These form factors are calculable on the lattice, but quite often prohibitively complicated if 
one uses staggered or Wilson fermions. The staggered fermions, with their weird flavor/taste struc- 
ture, make even definition of nucleon impractical. The Wilson fermions make the necessary cur- 
rent renormalization complicated due to explicit violation of chiral symmetry. The domain-wall 
fermions (DWF) [Q, with their exponential suppression of chiral symmetry breaking, make 

such renormalizations easy. In particular the vector and axial currents should share the same renor- 
malization. Thus the ratio g A /g v is naturally renormalized in DWF lattice calculations 

Structure functions are known from lepton deep inelastic scattering off nucleon the cross 
section of which is factorized in terms of the leptonic and hadronic tensors. The hadronic tensor, 
W^i V , is decomposed into symmetric unpolarized and antisymmetric polarized parts: 

w^(x,q 2 ) = ie" v "% (^(s,(x,e 2 )+s2(v,e 2 ))-^s2^e 2 )) 

with v = q-P, S 2 = -M 2 , x = Q 2 /2v. The unpolarized structure functions are F\{x,Q 2 ), F 2 (x,Q 2 ), 
and the polarized are g\ (x, Q 2 ), g 2 (x, Q 2 )- Their moments are described in terms of Wilson's oper- 
ator product expansion: 

2 f dx^F^Q 2 ) = £ c[% 2 /Q 2 ,g(Ll))(x n ) cl (Li) + 0(l/Q 2 ), 

Jo q=u,d 

[ l dxx"- 2 F 2 (x,Q 2 ) = £ c®fa 2 /Q?,g( f l))<f)M + 0(l/Q?), 
2 f 1 dxx n gl (x,Q 2 ) = £ e { fl^ 2 /Q 2 ,g(^)(x") Aq (^ + ff(\/Q 2 ), 

J0 q=u,d 

2 f dxx"g 2 (x,Q 2 ) = £ [^(AtVe^OOKOO 

q=u,d 

-24„( At 2 /G 2 , !? ( M ))(x") A ,( M )] + ^(i/e 2 ), 

where c\, c 2 , e\, and e 2 are the perturbatively known Wilson coefficients and {x n ) q {pL), (x"}^^) 
and d n {}i,) are calculable on the lattice as forward nucleon matrix elements of certain local oper- 
ators. Again the conventional staggered or Wilson fermions complicate such lattice calculations 
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for the same respective reasons as discussed about the form factors. The DWF calculations are 
simpler because of easier renormalizations due to good chiral symmetry. In particular the first mo- 
ments {x) u -d (quark momentum fraction) and (x)A H -Arf (quark helicity fraction) share a common 
renormalization and so their ratio is naturally renormalized in DWF calculations 

In this report we discuss dynamical DWF lattice QCD calculations of nucleon electroweak 
form factor ratio, g A /g v , the structure function ratio, {x} u -d/{x)Au-Ad, an d a polarized structure 
function moment d\ . The former two are naturally renormalized but d\ is yet to be renormalized. 



2. RBC and RBC/UK dynamical DWF ensembles 

The numerical calculations are performed with two different sets with dynamical domain-wall 
fermions (DWF) quarks. One employs two degenerate flavors of DWF quarks. The other employs 
three flavors, one strange and two degenerate and lighter up and down quarks. 

In generating the former set we used the rectangular-improved DBW2 gauge action [^, [1C| ]. 
The gauge coupling is set at 0.8 and the lattice cut off turned out to be about 1.7 GeV [O]. The 
four-dimensional lattice size is 16 3 x 32, corresponding to about (1.9fm) 3 spatial box. This set 
consists of three ensembles, each with different quark mass, 0.04, 0.03 and 0.02 in lattice units. 
These values roughly correspond to 1, 3/4 and 1/2 of physical strange quark mass. The pion mass 
for these dynamical quark mass values are about 700, 610 and 490 MeV, respectively. The nucleon 
mass are about 1.5, 1.4 and 1.3 GeV. The fifth-dimensional extent of the lattice is set as L s = 12 
and the domain-wall height M5 = 1.8. The residual mass was measured as m m& = 0.00137(5) or 
about 2.5 MeV. In this report we use 220 gauge configurations at the heavier two dynamical mass 
values and 175 at the lightest. These ensembles were generated using the QCDSP computers at 
RIKEN-BNL Research Center and Columbia University. 

The latter set was generated using another rectangular gauge action, the Iwasaki action [|f2|] . 



The gauge coupling is set at 2.13 and the lattice cutoff is about 1.6 GeV [ ]13| , 14J . There are two 
different four dimensional lattice sizes, 16 3 x 32 and 24 3 x 64, that corresponds to about (2.0fm) 3 
and (3.0fm) 3 spatial box respectively. In this report we concentrate on the latter, larger volume. 
The strange quark mass is set at about the physical value, 0.04 in lattice units. Three ensembles are 
generated with degenerate up and down quark mass set at 0.03, 0.02 and 0.01 respectively. The pion 
mass for these dynamical quark mass values are about 620, 520 and 390 MeV, the nucleon mass 
about 1.4, 1.3 and 1.2 GeV. The fifth-dimensional extent of the lattice is L s = 16 and the domain- 
wall height M 5 = 1.8. The residual mass was measured as m res = 0.00308(3) or about 4.8 MeV. In 
this report we use 25-30 gauge configurations at each up/down mass value. These ensembles were 
generated using the QCDOC computers at RIKEN-BNL Research Center, Columbia University 
and Edinburgh University. 



3. Lattice nucleon matrix elements 



Our lattice formulation follows the standard. The nucleon two-point function is defined as 
G N (?) = Tr[(l + Yt)Lx{TBi(x)Bi(0))} using B { = e ahc {u T a Cy 5 d b )u c for proton. The form-factor 
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three-point functions are defined as 

GY\tA = ^ + Yt)l^{TB l (xyr\x)B l {0))] 



and 



G«/(t,t>) 



1 



£ Tr[(l + 7 ,)W 5 II(rBi(x'Mr 4 WBi(0))] 

i=x,v,z x! x 



with fixed = f source — ? s ink set at about 1.5 fm in physical unit and varying t smaller than t' . The 

is j , lattice 

G»(t,t')-G d (t,t') 



ratio of their isovector combinations, 



directly yields the renormalized 



G-(t,t')-G d v {t,t')_ 
isovector charge ratio (g A /g v ) Ten - 

Unpolarized structure funtion moments are obtained from three-point functions, 



1 



EM'W-A.lM = 2 <0*(M)[Vtt. • • -P*, + • • • - (trace)] 



n-\ 



Ym Dpa'-'Dfr -(trace) 
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On the lattice we can measure (x) q , {x 2 ) q and (x 3 ) q . Higher moment operators mix with lower 
dimensional ones: operators belonging in irreducible representations of 0(4) transform reducibly 
under the lattice Hyper-cubic group. Only (x) q can be measured with P = 0. 
Polarized structure function moments are obtained by 

-Mtfow-H*}™ = ;^<*"MA0[^« 2 • • -P*, + ' ' ' " (^es)] 



2 ) T/HYoDh ■ D lln -(traces) 



(P^X^J^ = -TJ d M[(SoP^ ~SMP, 2 ■ ■ -P,„ + • • • - (traces)] 



f[5]q 



•x Y5Y[aD^] ■ ■ ■ Dn n -(traces) 



On the lattice we can measure (1)a? (&<0> {x)hq, (x 2 )A q , d\, and d^. Only (I)a^, {x)&q, and d\ can 
be measured with P = 0. 

The lattice bare values of these quantities need be renormalized before comparison with exper- 
iments. Such renormalizations are being performed non-perturbatively within Rome/Southampton 
RI/MOM scheme [|l5|, 16] and will be reported soon. In this report the two naturally renormalized 
ratios, g A /g v and (x) u -d/(x)Au-Ad, and unrenormalized d\ are discussed. 



4. Numerical results 

Let us first discuss the results obtained from the RBC two-flavor dynamical DWF ensembles. 
In figure |I] we plot the renormalized ratios, g A /g v (left) and (x) H _j / (x)au-M (right) against the bare 
up/down quark mass in lattice unit. Both show only mild quark mass dependence. g A /g v appears 
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Figure 1: Two-flavor dynamical DWF results (red cross) for renormalized ratios, g A /g v (left) and 
( x )u-d/ (x) Au-Ad (right) plotted against the bare quark mass. Respective experimental values are given as 
blue burst at m„ c / = — m res = —0.00137. Linear fit to g A /g v is also presented. 



in good agreement with the experiment. Linear fit to calculated values yields a value of 1.27(5) 
in the chiral limit, in comparison with the experimental value of 1.2695(29). (x) u -d / ' (x) Au-Ad also 
appears in broad agreement with the experiment. 

Now we turn to the results obtained from the RBC/UKQCD three-flavor dynamical DWF 
ensembles. These results are still preliminary with statistics of 25-30 configurations. In figure || 
we plot the renormalized ratios, g A /g v (left) and (x) u - c i / (x) Au-Ad (right) against the bare up/down 
quark mass in lattice unit. Again both show only mild quark mass dependence and appear in broad 
agreement with respective experiments. Linear fit to calculated g A /g v yields a value of 1.32(11). 

Since these calculated ratios are naturally renormalized, they can be compared with each other 
(see figure |3]) in physical scale such as pion mass squared, nr n . They appear mutually consistent 
and appear in broad agreement with the experiments. 

In figure Q we present the up and down contributions to the polarized structure function 
moment d\, with two (left) and three (right) dynamical flavors. This quantity summarizes the 
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Figure 2: Three-flavor dynamical DWF results (red cross) for renormalized ratios, g A /g v (left) and 
( x )u-d/ ( x ) Au-Ad (right) plotted against the bare up/down quark mass. Respective experimental values are 
given as blue burst at m uc / = — m res = —0.00308. Linear fit to g A /g v is also presented. 



5 



Nucleon structure with dynamical DWF 



Shigemi Ohta 



1.6 
1.5 
1.4 
1.3 
1.2 
1.1 
1 




-0.1 



0.1 



0.2 0.3 



0.4 



0.5 0.6 



1 

0.95 

0.9 
0.85 

0.8 
0.75 

0.7 
0.65 

0.6 
0.55 

0.5 




-0.1 



0.1 



0.2 0.3 0.4 0.5 0.6 



Figure 3: Renormalized ratios, g A /g v (left) and (x) u ^ l j/( x )Au-Ad (right) plotted against m\ in GeV 2 . Two- 
(+) and three-flavor ( x ) calculations and experiment (□). 



twist-3 part of the g2 polarized structure function. Wandzura-Wilzcek relation, g2{x) = —gi(x) + 
— gi(y), dictates it should be small [[17]]. However it need not be small in confining theories 



[18]. Our result, though yet to be renormalized, suggests it is small. 



5. Conclusions 



RBC two-flavor DBW2+DWF dynamical calculations are almost complete: the lattice cutoff 
is about 1.7 GeV, and the residual mass is m res = 0.00137(5) in lattice units. With three ensembles 
at degenerate up/down quark mass of m sea = 0.04, 0.03, and 0.02, the renormalized ratios, g A /g v 
and (x) u -d I (x) Au-Ad , appear in agreement with experiment, despite rather small volume. They 
show only mild quark mass dependence, and the linear chiral extrapolation for the former yields a 
value g A /g v = 1.27(5). Non-perturbative renormalizations for individual quantities are on the way. 
In particular, a polarized structure function moment, d\ , appears small, though not renormalized. 

RBC/UKQCD three-flavor, Iwasaki+DWF dynamical calculations are ongoing: the lattice cut- 
off is about 1.6 GeV and the residual mass is m res = 0.00308(3). With the strange quark mass of 
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Figure 4: Up (+) and down ( x ) contributions to bare polarized structure function moment d\ , with two 
(left) and three (right) dynamical flavors, against the bare up/down quark mass. They appear small. 
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0.04, three ensembles are being generated at up/down quark mass of 0.03, 0.02, and 0.01. Results 
from the large-volume (3-fm across) ensembles, though preliminary at small statistics, are encour- 
aging. The renormalized ratios, g A /g v and {x) u -d/ {x)& u -m, consistent with experiment. Again 
they show only mild quark mass dependence, and linear fit yields a value g A /g v = 1.32(11). d\ 
appears small again, though not renormalized. 
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were funded by the U.S. DOE grant DE-FG02-92ER40699, PPARC JIF grant PPA/J/S/1998/00756 
and by RIKEN. 

References 

[1] Yao, W.-M. et al. Review of Particle Physics. Journal of Physics G 33, 1+ (2006). 

[2] Kaplan, D. B. Phys. Lett. B288, 342-347 (1992). 

[3] Shamir, Y. Nucl. Phys. B406, 90-106 (1993). 

[4] Furman, V. and Shamir, Y. Nucl. Phys. B439, 54-78 (1995). 

[5] Blum, T. et al. Phys. Rev. D69, 074502 (2004). 

[6] Sasaki, S., Orginos, K, Ohta, S., and Blum, T. Phys. Rev. D68, 054509 (2003). 

[7] Taylor, R. E. Rev. Mod. Phys. 63, 573-595 (1991). 

[8] Orginos, K, Blum, T., and Ohta, S. Phys. Rev. D73, 094503 (2006). 

[9] Takaishi, T. Phys. Rev. D54, 1050-1053 (1996). 

[10] de Forcrand, P. et al. Nucl. Phys. B577, 263-278 (2000). 

[11] Aoki, Y. et al. Phys. Rev. D72, 1 14505 (2005). 

[12] Iwasaki, Y. unpublished UTHEF '-117 (1983). 

[13] Lin, M. PoS(LAT2006), 185 (2006). 

[14] Tweedie, R. PoS(LAT2006), 096 (2006). 

[15] Martinelli, G., Pittori, C, Sachrajda, C. T., Testa, M., and Vladikas, A. Nucl. Phys. B445, 81-108 
(1995). 

[16] Blum, T. et al. Phys. Rev. D66, 014504 (2002). 

[17] Wandzura, S. and Wilczek, F. Phys. Lett. B72, 195 (1977). 

[18] Jaffe, R. L. and Ji, X.-D. Phys. Rev. D43, 724-732 (1991). 



7 



